0-7r phase shifts in Josephson junctions as a signature for the s±-wave pairing state 
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We investigate Josephson junctions with superconducting ferropnictides, both in the diffusive and ballistic limit. 
We focus on the proposed s±-wave state, and find that the relative phase shift intrinsic to the s±-wave state 
may provide 0-7T oscillations in the Josephson current. This feature can be used to discriminate this pairing 
state from the conventional s-wave symmetry. The 0-7T oscillations appear both as a function of the ratio of the 
interface resistances for each band and, more importantly, as a function of temperature, which greatly aids in 
their detection. 
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The discovery of high-T c superconductivity in the 
ferropnictides^ has triggered an avalanche of investigations 
(see the reviews 2 and references therein) from a broad range 
of communities in condensed matter physics. A crucial issue 
which remains unresolved is the nature of the superconduct- 
ing order parameter (OP) symmetry in ferropnictide supercon- 
ductors. This topic is particularly intriguing since the ferrop- 
nictides feature a multiband Fermi-surface where the Cooper- 
pairs may reside. 

In order to identify the symmetry of the superconducting 
order parameter (OP), several recent experimental studies 3 4 
utilized the method of point-contact spectroscopy in order to 
study the symmetry of the superconducting OP in the ferrop- 
nictides. The findings were, however, not easily reconcil- 
able. Using an extended Blonder- Tinkham-Klapwijk (BTK) 
theory^ to fit their data, some groups 3 found a zero-bias con- 
ductance peak, indicative of a nodal order parameter such 
as d-wave. However, other groups 4 interpreted their data in 
terms of one or more nodeless OPs, such as s-wave. 

One of the leading candidates for the pairing symmetry is 
the so-called s±-wave state proposed in RefsP5 This pair- 
ing symmetry consists of two s-wave order parameters for the 
electron-like and hole-like Fermi surfaces that differ in sign. 
Some progress has been made in mapping out the ramifica- 
tions of the s±-wave sy mme try to quantum transport prop- 
erties of the ferropnictidesfflini p r instance, it has been pre- 
dicted that subgap surface states should appear in the presence 
of interband scattering 10 . Unfortunately, such subgap surface 
states are not unique for the s±-wave state, and do not provide 
unambiguous evidence for this pairing symmetry. 

To shed more light on the pairing symmetry in the ferrop- 
nictide superconductors, we present results for both the prox- 
imity effect and the Josephson current in hybrid structures in- 
volving normal metal elements and superconducting ferrop- 
nictides. The motivation for this is that both of these phenom- 
ena are expected to produce valuable information about the 
pairing state in the superconductor. We take into account the 
intrinsic multiband nature of this material class and include re- 
sults for the diffusive limit of transport, in contrast to previous 
theoretical works on these systems. 

For Josephson junctions with conventional superconduc- 
tors (s-wave), it is well-known that the supercurrent decays 
in a monotonous fashion as a function of both temperature 
and interlayer width, when the material separating the super- 



conductors is non-magnetic. If the interlayer is ferromag- 
netic, the current oscillates and goes to zero at certain critical 
widths and temperatures. This phenomenon is known as 0-7r 
oscillations 11 , and serves as a signature of either ferromag- 
netic correlations or nodal OPs, such as d-wave, present in the 
Josephson junction. 

In this Rapid Communication, we show that the aforemen- 
tioned prerequisites for 0-ir oscillations are rendered unneces- 
sary in the presence of an s±-wave pairing state. We find that 
0-7T oscillations may occur in a Josephson junction consist- 
ing of a conventional s-wave superconductor and a s±-wave 
superconductor separated by a normal (non-magnetic) inter- 
layer, thus in the complete absence of any ferromagnetic ele- 
ments or nodal superconducting OPs. This effect is explained 
in terms of the relative phase shift between the bands in the 
s±-wave superconductor, and constitutes a signature of the 
s±-wave state which can be probed in experiments. In fact, 
using such an observation in conjunction with other experi- 
ments that report a nodeless OP, ruling out d-wave pairing, 
would strongly support the realization of a s±-wave state. Our 
results are qualitatively independent of the interband scatter- 
ing strength, and are induced solely by the s±-wave symme- 
try. This renders our prediction more robust than recent pro- 
posals regarding subgap bound states as probes for the s±- 
wave state, which rely heavily on substantial interband scat- 
tering. 

We will employ the quasiclassical theory of superconduc- 
tivity in form of the Usadel 12 equation and the accompanying 
Kupriyanov-Lukichev boundary conditions^ modified for a 
multiband situational The quasiclassical approach is justi- 
fied under the condition that the Fermi energy is much larger 
than the superconducting gap and the impurity scattering self- 
energy, which should be a safe assumption for the ferropnic- 
tides. The notation and conventions of Ref. 15 will be used 
in what follows. For equilibrium situations, it suffices to con- 
sider the retarded part of the matrix Green's function, <?, which 
is parameterized conveniently by the quantity 9^, a =T,|- 
The Green's function satisfies g 2 = 1 and consists of entries 
with = cosh(d^) and = smh(d^) as measures of 
the proximity effect induced by the multiband superconduc- 
tor. In this parameterization, the Usadel equation's obtained 
as D^d^O 1 ^ + 2ies^ = 0, where Dn is the diffusion coeffi- 
cient in the normal metal and e is the quasiparticle energy. In 
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FIG. 1 : (Color online) Plot of the density of states (DOS) at x = (at 
the N|I interface) for a thin (djv/Cs = 0.2) and thick (d N /£s = 1.0) 
normal metal region. We have set ta = 0.5 and considered several 
values of r„. 




FIG. 2: (Color online) (a) Plot of the critical current for an s- 
wave|N| s-wave and s-wave|N|s±-wave junction, using r& = 1.0 
and |A S /Ai| = 1.0. (b) Plot of the critical current in the s- 
wave|N|s±-wave case, using |A S /Ai| = 0.5. In both (a) and (b), 
we have set djv/^s = 1.0. 



th e sup erconducting region, we use the bulk Green's functions 
g^ n l 15 lfor each band as denoted by the index A = (1,2), with 
belonging gaps A a = lA^le 1 ^. The unique feature of the 
s±-wave state is that the relative phase between the bands is 
7r, i.e. <pi = ip and if2 = f + tt, where <p is the superconduct- 
ing phase associated with the broken U(l) gauge symmetry. 

The Usadel equation must be supplemented with boundary 
conditions at the interface of the superconducting region. Un- 
der the assumption of a low interface transparency, we may 
employ generalized Kupriyanov-Lukichev boundary condi- 
tions which for an N|s±-wave interface at x = take the 
form d N g N d x g N \ x=dN = J2x ^[9n, 9\}\x=d N where d N is 
the thickness of the normal metal layer while 7^ = R^/Rn- 
Here, Rn is the resistance of the normal metal region, while 
Rg is the effective barrier resistance for band A. At x = 0, 
we have d x Q^ = 0, corresponding to zero outgoing current at 
the insulating/vacuum interface. 

Let us first briefly investigate the full proximity effect 
regime in a N|s± junction by solving the Usadel equation 
numerically with its boundary conditions. The normalized 
density of states (DOS) reads N(e)/N = E<x Re { c a ) 12 - 
There are three parameters that are free to vary in our the- 
ory. One is the thickness of the normal metal layer g?jv/£s, 
where £5 = y/D^ /\A\\. The two others are the ratio be- 
tween the gaps and the ratio between the barrier parameters, 
defined respectively as r& = IA2/A1J and r 7 = 72/71- In 
Fig. [T] we contrast the thin junction case d^/^s "C 1 with a 
thick junction djy/Cs — 1 for a representative choice of pa- 
rameters. We fix ta = 0.5 and plot the DOS in the N region 
at x = for several values of r 7 , with 71 = 5 correspond- 
ing to a low barrier transparency. There are in general three 
peaks in the energy-resolved DOS. Two of these peaks pertain 
to the bulk gaps of the s± superconductor, while the third de- 
marcates the opening of a minigap in the spectrum. This is 
qualitatively the same as what would be expected for a multi- 
band superconductor with a conventional s-wave symmetry, 
such as MgB;P. 

Therefore, the proximity effect and its impact on the DOS 
does not appear to provide a unique diagnostic tool in order 
to distinguish s±-wave symmetry from ordinary s-wave sym- 
metry. We thus turn our attention to the Josephson coupling 
for s±-wave superconductors as a possible mean to reveal this 
symmetry. To this end, we will consider a s-wave|N|s±-wave 



junction, where the s-wave gap is given by A s = \A s \e ltPs , 
and assume a weak proximity effect which allows us to lin- 
earize the Usadel equation and proceed analytically, facili- 
tating the interpretation of the obtained results. Also, the 
linearized approach is expected to yield excellent results in 
the experimentally relevant low-transparency case. The su- 
percurrent is given by Ij ~ J_ AeTx{py,{gd x g) K }, where 
f>3 = diag(l, 1, — 1, —1) and 'K' denotes the Keldysh compo- 
nent of the Green's function 11 . After solving the Usadel equa- 
tion, one may insert g into the above equation for the super- 
current. We find the following expression for the normalized 
zero-temperature Josephson current: 

poo 

I, j = Iq sin Aip, I = / dsRe{1l£/[ikdsm(kd)}}, (1) 
Jo 

where C = Ea^aM, K = Ea^J?M- Here ' 
Aip — (p — ip s is defined as the phase difference between 
band A = 1 in the right superconductor and the left super- 
conductor, k = ^J2ie/ Dn, while T^' R describe the anoma- 
lous Green's functions on the left/right side of the junction. 
These are proportional to the off-diagonal entries in the bulk 
Green's functions for the superconductors, which have the 
form !F\' R oc s^' R . We defined <5a=i = 1 an d <$a=2 = — 1- 
Note that the above expressions are valid for both a s-wave 
and s±-wave superconductor on either side of the diffusive 
normal metal, which is why we have included the band in- 
dex also on the left side. In the s-wave case, we have T\ = 
^sinh[arctanh(|A s |/e)], while in the s±-wave case we have 
T\ = sinh[arctanh(|A A |/e)]. 

We now solve Eq. ([T]l numerically to obtain the Joseph- 
son critical current, corresponding to I c = \Iq\ which is the 
relevant quantity measured experimentally. In Fig. |2|a), we 
plot the critical current as a function of the ratio between the 
interface barriers for each band, r 7 , for both a s-wave|N|s- 
wave and s-wave|N|s±-wave junction. In the former case, 
the current decays monotonously as is well-known. However, 
the situation is very different when we replace, say, the right 
s-wave superconductor with an s±-wave state. The current 
now displays 0-n oscillations, even in the complete absence 
of any ferromagnetic elements. This is very different from 
the conventional s-wave case, where a ferromagnetic element 
is required in order to induce the 0-ir oscillations. Thus, ex- 
perimental observation of such 0-tt oscillations in a Josephson 
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junction with ferropnictides would provide a strong indication 
of the presence of an s±-wave state. In Fig. |2|b), we give re- 
sults up to large r 7 for the s-wave|N|s±-wave case. As seen, 
the current saturates after the 0-ir oscillation since r 7 ^> 1 
means that one of the band interface transparencies tends to 
zero and does not contribute to transport. 

The appearance of the 0-7r oscillations in the current may 
be understood as follows. The transport of charge in an s- 
wave|N|s± -wave junction takes place both through inter- and 
intraband channels, as may be inferred directly by observing 
that the product C1Z in Eq. ([T]i produces precisely such terms. 
Due to the relative phase shift of tt between the two bands 
in the s±-wave state, these contributions to the critical cur- 
rent have opposite signs. For simplicity, consider the case 
where all gap magnitudes are equal in the Josephson junc- 
tion, |Aa| = |A S |, which leads to equal anomalous Green's 
functions T on both sides of the junction. We then have 
C1Z = T 2 (\j^\ — I/7D in Eq. {u, which is clearly seen 
to change sign at r 7 = 1. This does not occur in a conven- 
tional s-wave superconductor, where there is no relative phase 
shift. The basic mechanism behind the 0-tt oscillations is thus 
that variations in the barrier parameters "fxfor the bands will 
lead to either a dominant contribution between bands with no 
phase shift relative each other or bands with order parameters 
that differ in sign. 

Let us also consider the ballistic limit, to show that the 
mechanism for the 0-tt oscillations persists in clean samples. 
The only other change in the physical system under consider- 
ation is that we replace the normal interlayer with a thin insu- 
lating barrier (I), which in the BTK approach introduces the 
dimensionless barrier strengths Z\. In this manner, we can 
parameterize the relative barrier resistance in an analogous 
manner as with r 7 in the diffusive case by introducing rz = 
ZijZ\. We construct and solve the full 4x4 Bogoliubov- 
de Gennes equation for the two-band system, where we for 
generality also include coupling between the two bands pa- 
rameterized by the interband coupling strength a. This yields 
in general four current-carrying Andreev Bound States (ABS) 
E^(Atp). The Josephson current for this s-wave|I|s±-wave 
Josephson junction is then found in the ordinary way from 16 
I j = 2e£i = i where E l denotes the four ABS and 

f(E) is the Fermi-Dirac distribution function. 

To present an explicit illustration of the mechanism of 0- 
7r -oscillations in a s± -system in the ballistic limit, we pro- 
ceed analytically for the special case of a = 0. Here, 
we have for simplicity assumed that \A\\ = \A S \ = 
\A\. This gives solutions for the ABS on the well- 

knowrPform Ef = ±\A\^J I - D l sin 2 (A<^/2) and E% = 

±\A\y/l-D 2 cos 2 (Ay>/2), with D x = 4/(4 + Z 2 X ). At 
T = 0, the above expression for the Josephson current yields 
in the tunneling limit Ij = I\ sin ip, with I\ = (D\ —D-^IojA 
and I = 2e|A|. It is obvious that for Z 2 < Z\ one will have 
1\ < 0, i.e. the system being in the tt state, as explained for the 
diffusive case. As shown in Fig. [3] the crossover point above 
which the A = 1 contribution dominates instead is rz = 1. 
Notice however that the current does not vanish entirely at 
the crossover point due to a second harmonic component in 
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FIG. 3: (Color online) Critical current for a ballistic s-wavejljs±- 
wave Josephson junction as a function of the relative barrier strength 
rz- Interband coupling is neglected, and we have set Z\ = 6, T = 0, 
and \A\\ = |A S |. Inset: current-phase-relation for selected values 
of rz, as indicated by the arrows in the main figure. 



the current-phase-relation (as shown in the inset of Fig. [3]) 
which dominates close to the transition point. This is demon- 
strated explicitly by taking the approximation to the next order 
in the limit Z 2 = Z\, which yields Ij = I2 sin(2Ay), with 
I 2 = —IqD\/\6. We note that this non-sinusoidality of the 
current-phase-relation was absent in the diffusive case since 
the linearized Usadel equation corresponds to a first order ap- 
proximation in the interface resistance. We also emphasize 
that in this treatment, interband coupling is not essential for 
the occurrence of the 0-7r -transition. However, we have ver- 
ified numerically that the results of Fig. [3] are qualitatively 
valid also for a > 0, so that the predicted experimental signa- 
ture should be equally distinct for strong interband coupling. 

From the analysis above, it is seen that the crucial ingredi- 
ent for the observation of the 0-ir oscillations is having dif- 
ferent barrier parameters for each band A, or alternatively dif- 
ferent probabilities for Cooper-pair tunnelling. As suggested 
in RefP, these probabilities may be artificially altered by se- 
lecting materials with appropriate Fermi surfaces. Different 
Fermi vector-mismatches would then lead to different tunnel- 
ing probabilities. In our case, the size of the Fermi surface of 
the diffusive normal metal region could be modified by dop- 
ing. Thus, whereas 0-tt oscillations in S|F|S junctions can 
be seen as a function of the width dp of the ferromagnetic 
layei^l, necessitating the fabrication of several samples with 
different widths, the present scenario requires fabrication of 
several samples with the doping level in the normal metal 
varying in a systematic way. We note that it was also ob- 
served in Ref. 9 , although in the context of a superconducting 
s-wave | s± -wave | s-wave trilayer, that a 7r-junction could be 
fabricated in a similar manner. 

Although the above procedure is in principle feasible, it 
is very challenging to quantitatively relate the Fermi-vector 
mismatch directly to the parameter r 7 . However, we find 
that the 0-tt oscillations also occur as a function of tem- 
perature in the diffusive limit, thus constituting an alterna- 
tive, and simpler, approach to the recipe sketched above for 
altering r T . Assuming a Bardeen-Cooper-Schrieffer (BCS) 
temperature-dependence for the gaps, with a critical tempera- 
ture T c \ = T c for the s±-wave superconductor and TV, for 
the s-wave superconductor 20 , we plot the results in Fig. H As 
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FIG. 4: (Color online) Plot of the critical current as a function of 
temperature for an s-wavejNjs±-wave junction, using djv/^s = 1.0 
and |A S /Ai| = 0.5. In (a), we have r& = 0.3 while in (b) r& = 
1.3. 

seen, 0-ir oscillations appear as a function of temperature for 
a wide range of interface parameters r 7 . For large values of 
Ta, a normal monotonous decay of the critical current is seen. 
Although the exact relation between r 7 and r& which renders 
possible the Q-tt oscillations is difficult to extract analytically 
from Eq. ([T]), the basic mechanism is nevertheless the same 
as the one explained previously. From Fig. [4] we see that the 
absence of 0-tt oscillations not necessarily rules out that s± 
state, whereas the presence of them rules out the s-wave state. 

Finally, we point out that very strong impurity interband 
scattering T would eventually suppress the critical tempera- 
ture for the s± ground stated The difference between the 
DOS on the hole and electron Fermi pockets would determine 
how fast the suppression rate increases with T as compared to 
e.g. a d-wave scenario. For intraband scattering, however, the 
s± state is protected by Anderson's theorem. In our model, we 



have incorporated interband scattering only near the interface. 
A further extension of the model considered here could be 
to incorporate magnetic correlations in the s± state and also 
investigate strong interband scattering in the bulk of the su- 
perconductor to see how it affects the transport properties^, 
although we expect that they would remain qualitatively the 
same as reported here since the basic mechanism for the 0-7r 
oscillations would remain intact. 

In summary, we have investigated the Josephson coupling 
properties of junctions with s±-wave superconductors. In 
contrast to previous literature, we have here included results 
for both the ballistic and diffusive regimes. The relative phase 
shift of the bands intrinsic for the s±-wave state leads to 0- 
7r oscillations in an s-wave|N|s±-wave Josephson junction, 
even in the absence of any ferromagnetic elements. The mech- 
anism behind these oscillations is a competition between the 
sign-dependent contribution of transport from different bands 
in the s±-wave superconductor to the s-wave superconductor. 
The 0-7r oscillations are seen as a function of temperature, thus 
vastly facilitating the experimental testing of our predictions 
compared to methods that involve changing the parameters 
of the model system. Our results may aid in identifying the 
possible existence of an s±-wave pairing state in the super- 
conducting ferropnictides. 
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